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Strange attractor simulated on a quantum computer 
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Abstract. We show that dissipative classical dynamics converging to a strange attractor can be simulated 
on a quantum computer. Such quantum computations allow to investigate efficiently the small scale struc- 
ture of strange attractors, yielding new information inaccessible to classical computers. This opens new 
possibilities for quantum simulations of various dissipative processes in nature. 

PACS. 05.45.Df Fractals - 05. 45. Ac Low-dimensional Chaos - 03.67.Lx Quantum Computation 



Starting from the work of Lorenz H, it has been re- 
alized that the dynamics of many various dissipative sys- 
tems converges to so-called strange attractors ||. These 
objects are characterized by fractal dimensions and chaotic 
unstable dynamics of individual trajectories (see e.g. ||[|). 
They appear in nature in very different contexts, includ- 
ing applications to turbulence and weather forecast ^,^| , 
molecular dynamics Q , synchronization || , chaotic chem- 
ical reactions , multimode solid state lasers j|] an d com- 
plex dynamics in ecological systems ||,[l0) and physiology 
[jllj . The efficient numerical simulation of such dissipative 
systems can therefore lead to many important practical 
applications. 

Recently, it has been understood that quantum me- 
chanics allows to perform computations in a fundamen- 
tally new way (see for review e.g. | p^ , p^l4]Jl5| , ^6t ) . Indeed, 
quantum parallelism can enormously accelerate the com- 
putation and provide new information inaccessible to clas- 
sical computers. Well-known examples are Shor's factor- 
ization algorithm p7) , which is exponentially faster than 
any known classical method, and Grover's search algo- 
rithm |l8| ], where the gain is polynomial. Even if impor- 
tant progress has been achieved during the last years, still 
it is essential to find new areas where quantum processors 
might give access to new information unreachable classi- 
cally. Especially interesting are applications to dissipative 
systems with irreversible dynamics leading to a loss of in- 
formation. 

In this Letter we analyze how classical dissipative dy- 
namics can be simulated on such quantum processors. To 
this aim, we study a simple deterministic model where 
dynamics converges to a strange attractor, and show that 
it can be efficiently simulated on a quantum computer. 
Even if the dynamics on the attractor is unstable, dis- 
sipative and irreversible, a realistic quantum computer 
|0,|13|Jl4|,na,n9 20 can simulate it in a reversible way, and, 



To study how a quantum computer can simulate dissi- 
pative dynamics leading to a strange attractor, we choose 
the deterministic map given by: 



y = y/2 + x (mod 2) 



y (mod 1) . (1) 



already with 70 qubits, will provide access to new infor- 
mations inaccessible for modern supercomputers. 



where —0.5 < x < 0.5, — 1 < y < 1 and bars note the new 
values of variables. The map has one positive A+ and one 
negative A_ Lyapunov exponents (A± = ln[(5 ± \/l7)/2]) 
so that the dynamics converges to a strange attractor 
with Hausdorff (Dh) and information (-D/) dimensions 
[|l]j2|,p D H « D[ = 1 + A+/|A-| w 1.543. 

To implement this map for a computer simulation it 
is necessary to discretize the phase space. We choose the 
natural discretization in the binary representation of co- 
ordinates (x, y), so that the dynamics takes place on a 
regular square lattice with TV x 2N points, with N = 2 Uq 
and n q integer. The division in ([l]) is realized by shift and 
truncation of the last binary digit. With this procedure, 
the dissipation generates a discretized irreversible map, 
displaying a discretized strange attractor (see Fig.l top) 
which approaches the continuous one for large N . Such a 
map can be implemented efficiently on a quantum com- 
puter. For that, an initial image with Nd points is coded in 
the wave function |<?o >= Y^i j a ij\ x i > \Uj > |0 > |0 >, 
where a^j — or 1/y/Nd- Here the two registers \x{ > 
and \yj > with n q and n q + 1 qubits hold the values of 
the coordinates x and y (xi = —0.5 + i/N, i = 0, N — 1 
and yj = -1 + j/N,j = 0, ...,2N - 1). The third register 
with n q — 1 qubits is used as workspace for modular ad- 
ditions, and the last one collects the truncated last digits 
generated by the divisions ( "garbage" ) . We start with the 
simplified algorithm for which the garbage gets one digit 
at each map iteration so that t iterations need t qubits 
in the fourth register. The size of this register can be sig- 
nificantly reduced using a more refined algorithm we will 
describe later. 

The algorithm starts with the initial state >; first 
it places the last qubit of \y > in the garbage register, and 
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Fig. 1. (color) Top: exact classical/quantum computation of 
the discretized map ([!]); initial image (left) converged to the 
strange attractor after t = 10 iterations (right). Middle: quan- 
tum computation with noise amplitude e = 0.05 in each gate 
operation; attractor at t — 10 (right) and initial image recov- 
ered after 10 backwards iterations with fidelity / = 0.63 (left). 
Bottom: same as middle with e = 0.1 and / = 0.15. Left shows 
the central cell (—0.5 < x,y < 0.5), right shows the whole 
phase space. Color marks the probability density (integrated 
over third and fourth registers), from blue (density less than 
10~ 5 ) to red (maximal value). Here n q = 6, with in total 28 
qubits used. 



uses n q swap gates to shift the \y > qubits and obtain 
y/2. Then a modular addition is implemented in the way 
described in 0] to add x to y/2. After that, another mod- 
ular addition adds the second register to the first. In total, 
this requires 17n q — 10 quantum operations using Toffoli, 
control-not and swap gates in contrast to 0(2 2 "«) opera- 
tions for the classical algorithm. It is interesting to note 
that the algorithm allows to restore the initial state: in- 
verse map iterations are performed (a; = x— y, y/2 — y—x) 
and y is restored from y/2 using a qubit stored in the 
garbage register. This requires a similar number of oper- 
ations as the forward iterations. In principle this can be 
done on a classical computer in 0(2 2rlq ) operations, but 
this requires additional exponentially large memory which 




Fig. 2. (color) Spectral density \C{t, k Xty )\ 2 for the strange 
attractor of Fig.l at t = 10 in the region -TV/2 < k x < N/2, 
— TV < k y < N (TV = 2™«). Top: three-dimensional plot of 
the full distribution for n q = 6 and e = (left), e = 0.025 
(right). Middle: color plot of the coarse-grained distribution 
with n q = 6, rif =4 and e = (left), e = 0.025 (right). 
Bottom: same as middle for n q — 10, e = (left), and n q = 6, 
e — 0.05 (right). Colors are as in Fig.l. 

stores about 2 2r M bits for TV^ ~ TV 2 , contrary to only t 
qubits used by the quantum computer. 

Fig.1 shows the dynamics generated by the discretized 
map ([0) simulated on a quantum computer with exact 
and noisy unitary gates with imprecisions of amplitude 
e. Due to the dissipative nature of the map (|l|), the ini- 
tial image rapidly converges towards the strange attractor 
(already t = 5 is enough). Even in the presence of rela- 
tively strong noise, the fractal structure of the attractor is 
well-preserved and the initial image can be reliably recov- 
ered after backwards iterations, despite the exponential 
instability of the classical dynamics (|l]). The precision of 
computation can be quantitatively characterized through 
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the fidelity / defined at a given moment of time as the 
projection of the quantum state in presence of gate im- 
perfections on the exact state without imperfections. The 
global properties of the initial image can be recovered even 
at relatively low fidelity values. 

Even if the quantum algorithm performs one map iter- 
ation only in 0(n q ) operations, it is important to take into 
account the measurement procedure that allows to extract 
efficiently the information coded in the wave function. In- 
deed the number of points in Fig.l grows exponentially 
with n q and an exponential number of measurements is 
required to obtain the full density distribution. However, 
certain characteristics can be extracted in a polynomial 
number of measurements, providing new information in- 
accessible for classical computation. An example of such a 
quantity is the spectrum of phase space correlation func- 
tions, defined as C(t, k x , y ) = ^2 Xo , yo exp(2in(x(t, x , y ) + 
y(t, xq, yo))) exp(2iir(k x xo + k y yo)), where the sum runs 
over the points (xo,yo) of the initial distribution, and 
(x(t,x ,y ),y(t,x ,y )) is the position of {x ,y ) at time 
t. Such correlation functions have been studied for chaotic 
systems, where they determine various kinetic coefficients, 
for example the diffusion rate (see e.g. Q] p. 328). Due 
to chaos, the function C(t,k x ^ y ) has significant values at 
exponentially high harmonics k XiV ~ exp(|A_|t) which 
rapidly reaches harmonics of order N. In the theory of 
classical chaotic dynamics it is well-known that the infor- 
mation about such harmonics is very hard to access, since 
exponentially small scales should be explored, which can 
be done only with exponentially many trajectories p|,p5||. 
On the contrary, the quantum computation of C(t,k x , y ) 
can be done efficiently. For that, one makes t iterations of 
(Q), and creates the state a x,y exp(2in(x + y))\x > \y > 
|0 > \g(t) >. The preparation of this state is easily done by 
applying 2n q + 1 one-qubit rotations to the first two regis- 
ters. Then the garbage g(t) is erased by iterating the map 
backwards t times, that at the same time returns the coef- 
ficients a x , y to their original values. This creates the state 
E a ni* exp(2i7r(x(£,a;o,?/o) + y(i,2;o,2/o))|iz:o > \y > |0 > 
|0 >, keeping phases unchanged. The whole procedure is 
sometimes called "phase kickback". After that the appli- 
cation of a two-dimensional quantum Fourier transform 
Ipj yields in 0(n 2 ) operations the state Yl C(t, k XlV )\k x > 
\k y > |0 > |0 >. A polynomial number of measurements 
yields the principal peaks, or enables to obtain a coarse- 
grained image of the spectral density \C(t,k Xttl )\ 2 in the 
Fourier space. Indeed, independently of n q , one can mea- 
sure the first rif and nf + 1 qubits of the x and y regis- 
ters respectively, that gives integrated probability inside 
2 2n /+ 1 cells. Fig. 2 displays the spectral density for the case 
of Fig.l. It shows that new information about the coarse- 
grained spectral density can be obtained efficiently. In- 
deed, patterns are clearly present in Fig. 2 and they vary 
irregularly with n q (compare Fig. 2 middle left and bot- 
tom left). This confirms the nontrivial nature of informa- 
tion provided by the coarse-grained density \C(t, k XtV )\ 2 . 
Although the spectral density is more sensitive to noise 
than the distribution in Fig.l still the patterns remain 
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Fig. 3. Complexity of classical and quantum computations 
of coarse-grained distribution of Fig. 2 with 10% accuracy (i.e. 
fidelity f(t) = 0.9). Here n/ = 4, t = 10 and different values 
of n q are shown. For classical computation (open circles), M 
gives the number of classical Monte Carlo trajectories required 
to reach 10% accuracy. For quantum case (full circles), M is the 
number of measurements needed to obtain the same accuracy. 
Full line shows the total number of classical points in the initial 
image shown in Fig.l (M ~ 0.36 x 2 2 ™ 9 ). Logarithm is decimal. 



well-defined even in the presence of relatively strong er- 
rors (see Fig. 2). 

It is important to stress that even modern supercom- 
puters are unable to find the properties of the spectral 
density for n q > 20. Indeed, as is shown in Fig. 3, a classi- 
cal Monte Carlo algorithm requires an exponentially large 
number of trajectories M (M = 0(2 2riq )) to obtain the 
coarse-grained spectral density at fixed nf with fixed ac- 
curacy. In contrast, the quantum computation requires a 
number of measurements M independent of n q (each mea- 
surement is done after t map iterations and one Fourier 
transform which needs 0{n 2 ) quantum gates). 

To study the effect of noisy gates in a more quantita- 
tive way, we show in Fig. 4 the dependence of the fidelity 
f{t) of the quantum computation of spectral amplitudes 
C(t,k XlV ) on the noise amplitude e and total number of 
gates n g applied [n g = t{AAn q - 14) + (n q + 2) 2 - 2). The 
data show the global scaling law 1 — f(t) « e 2 n g /lA, valid 
for moderate e. The physical origin of this scaling law is 
related to the fact that for randomly fluctuating unitary 
gates the loss of probability from the exact state is of the 
order of e 2 for each gate operation. This law determines a 
time scale tf ~ l/(6e 2 n q ) up to which a reliable quantum 
computation is possible (/ > 0.5). Beyond tf the decoher- 
ence destroys the accuracy of the quantum computation 
and the results become strongly distorted, see e.g. Fig. 2 
bottom right. We note that a similar time scale appears 
in quantum computation of Shor's algorithm on a realis- 
tic quantum computer J2(| . The computation beyond the 
scale tf is possible but requires the application of quan- 
tum error-correctin g co des, at the cost of additional qubits 
and gates (see e.g. [ |14jjl5|l and refs. therein). Without er- 
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Fig. 4. Fidelity / for the spectral amplitudes of Fig. 2 as a 



function of e"n„ 



< 6, 10" 4 < e < 0.5 and 6 < t < 10. 



Inset: error probability W g in the third and fourth registers as 
a function of e 2 t for n q = 6. Straight lines show the theoretical 
slope 1, logarithms are decimal. 



ror correction, it is still possible to improve the fidelity of 
the final state by measuring the third and fourth regis- 
ters. Indeed for exact computation they are at zero after 
the backwards iterations while with noisy gates the er- 
ror probability grows as W g « 3.5e 2 t (Fig. 4 inset). The 
measurements of these registers allow to select the correct 
states and increase the fidelity by a factor 1/(1 — W g ), e.g. 
for the case of Fig.l bottom this procedure gives / = 0.22. 

The above algorithm is optimal for not very large times 
t. If one is interested in simulating the dynamics on the 
attractor for large t, then the size of the garbage regis- 
ter can be significantly reduced. Indeed, at any t, copying 
the results x{t),y{t) in two additional registers and re- 
versing the sequence of gates allows to erase the garbage 
and reproduce Xo,yo. This procedure can be done recur- 
sively following the strategy of "reversible pebble game" , 
the description of which can be found in Jl5[ . In the sim- 
plest version, n t qubits in the garbage register (plus the 
two additional registers for Xo,yo) allow to perform map 
iterations up to t ~ 2™* . This gives only a polynomial in- 
crease in the number of elementary quantum operations, 
being proportional to n^ 58 . The procedure becomes cost- 
effective for t ^> n q . 

The algorithm described above can be generalized to 
other dissipative maps. For example, modular multiplica- 
tions can be performed in 0(n 2 ) operations, as described 
in [E4| . This allows to simulate efficiently the map (0) with 
x 2 term in first equation and also the Henon attractor 
x — y + 1 — ax 2 , y — bx j4|] . Such an algorithm can also 
be adapted to perform the finite-step integration of the 
Lorenz system x = —a(x — y), y — ~xz + rx — y, z = 
xy — bz Jy . However the simulation of the dissipative dy- 
namics of these models requires more qubits than for (|l|) . 

Thus on the example of the map (Q), we have shown 
that a quantum computer can efficiently simulate dissi- 
pative irreversible dynamics. A quantum processor with 



70 qubits will be able to provide new information about 
small-scale structures of strange attractors inaccessible to 
modern supercomputers. 
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